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Abstract — In this study, towards analog and digital 
elliptic filters, use of transfer function through the choice 
of a particular filter structure from a set of design 
specifications to achieve the step change in transfer 
function analysis. Elliptic filter design with two 
optimization criteria (1) maximally flat amplitude 
response and (2) small amplitude ripple response, and 
the excellent performance of the filter available. There 
are many different filter design methods discussed in this 
paper is the first for the analog method, used in analog 
filters. The specified frequency response with amplitude 
response and phase response which obtained through the 
rational transfer function approximation. In addition to 
the order of the filter for different analysis, the use of 
continuous signals into discrete Z-transform system by 
using (1) Pulse constant conversion method (2) Bilinear 
transformation method two conversion methods. This 
paper and using MATLAB software to calculate the 
analysis results, and the results plotted graphically 
confirmed that the two algorithms are calculated from 
the Z-transform is very similar to results with the 
correct solution, using MATLAB software is used to 
evaluate the convenience, and By the results of its 
analysis to support the correctness of the conversion 
result, and confirmed that the bilinear transformation 
method than the pulse of the same conversion method 
has a better conversion results. 

Keywords — Audio and Video, Elliptic filter, Signal 
processing, Matlab, Bilinear transformation method, 
Impulse invariance method 

I. Introduction 

Elliptic filter widely used in the realization 
of speech processing. Because in speech 
analysis and speech synthesis that has the 
excellent processing power ability, the 
mathematical model of analysis by a transfer 
function can be described, or you can use the 
amplitude response, phase response 
described. Through literature Ma et al. [1] 
applied to the circuit signal analysis, Romero 
et al. [2] engaged in the ladder filter of the 
test, Gordon et al. [3] the reduction in filter 
design and analysis of their response, 
Alarcon et al. [4] the use of three-pole 
Butterworth analog filter for ECG signal of 
the research. The study of algorithms [5], 
first to use the Laplace transform, from the 
regularization of the classical analog low- 
pass filter. For the pre-distorting the cutoff 



frequency, so the other frequency selective 
filter is also applicable, such as high-pass, 
band pass and band stop filter. For the 
frequency conversion of analog design with a 
diameter cut-off frequency at different 
degrees of order low-pass Butterworth filter 
is normalized Butterworth low-pass filter 
start and uses the filter cutoff frequency and 
order of replacement which is s functional. 
This is an example of frequency conversion. 
Using the same method can be normalized 
low-pass filter, converted to other types of 
frequency selective filters. Finally, the output 
of the system to produce a continuous time 
signal reconstruction filter, on behalf of the 
original input signal is filtered results. This 
paper and using MATLAB software to 
analyze the results of calculations [6], and 
the results graphically plotted for reference 
comparison. Since elliptical filters better than 
the Butterworth filter's performance, so this 
article gives special seminars. Recently a 
more linear elliptic filter, designed to be 
more flexible and has a wide bandwidth [7,8], 
on the other filter current transmission 
control is also a new development [9], and 
have very low band-pass sensitivity and 
compensation effect [10], usually require 
high-order elliptical filter before they can 
achieve the above effect [11-14]. Through 
references [1-14] and to consolidate research 
and discussion to make a design of high- 
order elliptic filter is one of the points this 
study, this study first use of the Laplace 
transform algorithm, followed by its digital 
technology, making it the digital filter and 
use the calculation (computation) to perform 
continuous-time signal filtering action. 
Design frequency-selective filter, is used to 
convert the continuous-time signal 
corresponding to the series. Digital filter in a 
series of sampling points on the basis of 
processing sequence to generate new series, 
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new series followed by a digital to analog (D 
/ A) converter, conversion to the 
corresponding continuous time signal. 
Finally, the output of the system to produce a 
continuous time signal reconstruction filter, 
on behalf of the original input signal is 
filtered results. In this study using Impulse 
invariance method and Bilinear 
transformation method to deal with the 
digital signal transfer, change of continuous 
signals into discrete Z-transform system. 

II. Main Contains 

Signal processing algorithms often occur in 
a series of complex operations can easily 
make the calculation error, re-verification 
must spend a lot of time in order to calculate 
the convenience and accuracy of using 
MATLAB to write program makes saving a 
lot of time on operations, and therefore Our 
main areas of this monograph is often used in 
engineering to math up development, such as; 
linear algebra, polynomial solving, 
probability calculations ... etc., developed the 
program features allow users more 
convenient operation. Elliptic filter is 
sometimes called Cauer filter. This filter has 
the same ripple for band pass and band stop. 
Discussed in the previous filter in the known 
order of the filter, band pass and stop with 
error, Elliptic filter with minimal transient 
bandwidth. The size of the square of the 
response is 

|H(jO)| = j—! 

(1) 

Where U N (x) is the N order elliptic function, 

while s is the pass band ripple of the 
constant. 

2-1 Frequency transformation 

There are four basic types of frequency 
selective filters, for low pass, high pass, band 
pass and band stop. Discussed in the previous 
design skills, we only take into account the 
low-pass filter. Low-pass filter can be 



regarded as the prototype filter and the 
system functions fairly easy to obtain. 
Furthermore, to design high-pass or band- 
pass or band stop filter, frequency conversion 
can be used easily to obtain. Frequency 
conversion can be done in two ways. In the 
analog frequency conversion, the analog 
prototype filter h p (s) into a desired system 

function of the digital filter of another analog 
system function H(s) . Then use the 
corresponding skills, into a system 
function H(z) of the digital filter. In the 
digital frequency conversion, the analog 
prototype filters into a first digital domain, 
with the system function h p (z)- Then use of 

frequency conversion, so that can be 
converted into the digital filter. 

2-2 Analog frequency transformation 

Analog frequency conversion formula used 
in a prototype low-pass filter into a low-pass 
(with a different cut-off frequency), high 
pass, band pass or band stop filter, as shown 
below. 

(i) low-pass filter cut-off frequency q, to a 
new low pass filters cut-off frequency^ *: 

Qc 

s ^7^ 

(2) 

Therefore, if the prototype filter response of 

the system isH p (s), the low -pass filter of the 

system response will 
H(s) = H P — 

(3) 

(ii) low-pass filter cut-off frequency f\, to a 

new high-pass filters cut-off frequency £\ * : 

s 

(4) 

the system function of high-pass filter is 
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H(s) = H p ^ 

I s J 

(5) 

(iii) low-pass filter cut-off frequency Q c to 

the cut-off frequency Q , and high cut-off 

frequency of the band pass filter Q 2 : 



• a 



s +fi ,Q 2 
s(Q !-Q 2 ) 



(6) 

the system function of band-pass filter is 



Qr 



s 2 +QiQ 2 A 



H(s) = H f 

(V) 

(iv) cutoff frequency Q for low pass filter 

cutoff frequency to the low Q , and high 

cutoff frequency Q 2 of the band stop filter: 

s(Q 2 - fi i) 



.Q 



s +Q ,Q 2 



(8) 

the system function of band-stop filter is 



H(s) = H, 

(9) 



Q 



s(Q 2 -Q,)' 
c s 2 +Q ,Q 2 



2-3 Digital frequency conversion 

As in the analog domain, frequency 
conversion in the digital domain is feasible. 
In the digital domain into the frequency 
converter by replacing the function Z" 1 to 
achieve, that is, f(Z~'). This corresponds to 
take into account the stability criteria. All 
poles must fall within the unit circle 
corresponds to its own, and also corresponds 
to the unit circle itself. To achieve the 
corresponding unit circle to itself, for r = 1 , Z 



r ""e 



e -J ra = f ( e - jm ) = f ( e - Jro ) 



j'arg[f(e- i ")] 



(10) 



Analog frequency transformation formulas as 
shown in Table (1) 



Therefore , to the all frequency , must 
be \f( e - jm )\ = \ o So that, the 

correspondence all pass filter can be written 



as 



f( Z - i )=+n 



k=ll-Ot k Z 



(11) 



Table (1) Analog frequency transformation formulas 



Analog Transformation 



T Qc 

Low pass s -^ — r s 



High pass s -^ 






Band 



f s 2 +Q,n, A 



Band 



s->Q 



s(Q 2 _ ^l) 

s(Q 2 - Q i) 

C ? 

s +Q jQ 2 



J 



pass 



stop 



In order to stabilize the prototype filter to be 
stable filter must set |ak| < 1 . For the 

conversion prototype low-pass digital filter 
to the digital low-pass, high-pass, band-pass, 
or band-stop filter, the conversion formula 
from (10) obtained, as shown in Table (2). 



Table (2) Conversion prototype low- 
pass DIGITAL FILTER TO THE DIGITAL LOW- 
PASS, HIGH-PASS, BAND-PASS, OR BAND-STOP 
FILTER 



Types 


Transformat 
ion 


Design 
parameters 


Low- 
pass 


, z" 1 -a 

z -> r 

1-az" 1 


s'mUm -co') 12] 

a = — — — 7- — - 

s'm[(a> c + co" c )l 2] 


High- 
pass 


, z -1 +a 

z -* r 

1+az" 1 


cosUa) -<y*)/2] 

a ^ 

cos[(oi c + & c ) 1 2] 
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Band- 


i z 2 -qz'+a 
z — > 


a, =-2aJC/(K+l) 


pass 


<^z~ 2 -qz~'+l 


a, =(K-\)l(K+\) 

cos[(cc^+ai|)/2] 
cos[(co c -cq)/2] 


Band- 


,-i ^-of'^h 


a, = -2a/(K+l) 


stop 


ap -(\z +1 


a 2 =(l-K)/(l+X) 

cos[(ft), +fti )/2] 

a= ! 

cos[(a>,-ffi|)/2] 



Magnitude Response 



Frequency conversion can be done either two 
available techniques, however, important to 
note that what techniques to use. Move out, 
the pulse invariant transformation is not 
suitable for high frequency high-pass 
resonant or band-pass filter. In this case, 
assuming a low-pass prototype filter which 
uses the analog frequency conversion into 
high-pass filter, and then using impulse- 
invariance techniques into the digital filter. 
This will cause aliasing problems. However, 
if the same low-pass filter prototype, the first 
use of pulse techniques to convert non- 
variant digital filter, then after using the 
digital frequency conversion into high-pass 
filter, phase will not have any aliasing 
problem. Regardless of when to use the 
bilinear transformation, whether using analog 
or digital frequency conversion frequency 
converter is not important. In this case, 
analog and digital frequency conversion 
techniques will get the same results. 



in. MATLAB Examples 

Example 1 : In accordance with the following 
specifications makes analog elliptical low- 
pass filter design. Wp = 0.27t; Ws = 0.3ti; 
Rp=ldB; As=16dB; 




Analog frequency in pi units 
Magnitude in dB 




Analog frequency in pi units 
Phase Response 



Analog frequency in pi units 
Impulse Response 




time in seconds 



Figure 1. The relationship for analog elliptical low-pass filter 
frequency to amplitude and phase 



From above Example 1 we can find out that 
analog elliptical low-pass filter at the confine 
specifications Wp = 0.2ti; Ws = 0.37t; Rp = 
ldB; As = 16dB will produce a small ripples 
at the low pass region and obtain a clearly 
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slope for the next region, it is a good 
transform results for a low pass elliptic filter. 

Example 2: In accordance with the following 
specifications makes digital elliptical low- 
pass filter design. Wp = [OAn O.671]; Ws = 
[0.3ti 0.75tc] ; Rp = ldB; As = 40dB; 
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Figure 2. The relationship for digital elliptical low-pass filter 
frequency to amplitude and phase 

From above Example 2 we can find out that 
analog elliptical band-pass filter at the 
confine specifications Wp = [0.4tc 0.6tc]; Ws 
= [0.3tc 0.75tc] ; Rp = ldB; As = 40dB will 
produce a small ripples at the band pass top 
region and obtain a clearly slope for the other 
regions, it is a also good transform results for 
a band pass elliptic filter. 

Example 3: In accordance with the following 
specifications makes digital elliptical low- 
pass filter design. Wp = 0.3; Ws = 0.4; Rp = 
ldB; As = 20dB; 

ellip filter frequency response 
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Figure 3-a. The relationship for digital elliptical low-pass filter 
frequency to amplitude and phase 
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Figure 3-b. Digital elliptical low-pass filter pole and zero points 
map 

From above Example 3 we can find out that 
digital elliptical low-pass filter at the confine 
specifications Wp = 0.3; Ws = 0.4; Rp = 
ldB; As = 20dB will produce a small ripples 
at the low pass region and obtain a clearly 
slope for the next region, it is a good 
transform results for a low pass elliptic filter. 
On the other hand, from Figure 3-b has 
shown that all the pole and zero points are 
convergence at the unit circle. 

Example 4: In accordance with the following 
specifications makes digital elliptical high- 
pass filter design. Wp = 0.4; Ws = 0.3 ; Rp = 
ldB; As = 20dB; 

ellip filter frequency response 
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Figure 4-a. The relationship for digital elliptical high-pass filter 
frequency to amplitude and phase 




Figure 4-b. Digital elliptical high-pass filter pole and zero points 
map 

From above Example 4 we can find out that 
digital elliptical high-pass filter at the confine 
specifications Wp = 0.4; Ws = 0.3 ; Rp = 
ldB; As = 20dB will produce a small ripples 
at the high pass region and obtain a clearly 
slope for the next region, it is a good 
transform results for a high pass elliptic filter. 
On the other hand, from Figure 4-b has 
shown that all the pole and zero points are 
convergence at the unit circle. 

Example 5: In accordance with the following 
specifications makes digital elliptical Band- 
pass filter design. Wp = [0.3 0.5]; Ws = [0.2 
0.55]; Rp=ldB; As = 20dB; 

ellip filter frequency response 




0.4 0.6 

freq (rad/sec) 



0.4 0.6 

freq (rad/sec) 



Figure 5-a. The relationship for digital 
elliptical band-pass filter frequency to 
amplitude and phase 
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Figure 5-b. Digital elliptical band-pass filter 
pole and zero points map 

From above Example 5 we can find out that 
digital elliptical band-pass filter at the 
confine specifications Wp = [0.3 0.5]; Ws = 
[0.2 0.55] ; Rp = ldB; As = 20dB will 
produce a small ripples at the band pass 
region and obtain a clearly slope for the next 
region, it is a good transform results for a 
band pass elliptic filter. On the other hand, 
from Figure 4-b was shown that all of the 
pole and zero points are convergence at the 
unit circle. 

Example 6: In accordance with the following 
specifications makes digital elliptical Band- 
stop filter design. Wp = [0.2 0.55]; Ws = 
[0.3 0.5] ; Rp = ldB; As = 20dB; 

ellip filter frequency response 
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Figure 6-a. The relationship for digital 
elliptical band-stop filter frequency to 
amplitude and phase 



Heai ran 

Figure 6-b. Digital elliptical band-stop filter 
pole and zero points map 

From above Example 6 we can find out that 
digital elliptical band-stop filter at the 
confine specifications produce a good 
transform results for elliptic filter. The most 
important thing is finding that the bilinear 
transformation method is better than the 
Impulse-invariance method for the high pass 
filter status. 

Example 7. Impulse-invariance method and 
bilinear transformation method will be 
converted analog elliptical low-pass filter to 
high-pass digital filter. Filter specification 



are 



(B pL = 0.271, fi) pH 



ldB 



0.2rr,R p 
co s =0.3tc,A s =15dB 
ra pL :Low-pass digital frequency 
co pH : High-pass digital frequency 
R p : Ripple of band-pass 
co s : Band-stop digital frequency 
A s : Band- stop digital frequency decay 

This article by the above theoretical analysis, 
first create a simulation system by Laplace 
transform method to expand as part of the 
fractions, and using MATLAB software to 
analyze the selected time interval of its 
calculation and the results plotted graphically 
in Figure 1 . Impulse-invariant method and 
Bilinear transformation method will be 
converted to analog elliptical low-pass filter 
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and high-pass digital filter. 



pass Filter Magnitude Respon: 




frequency in pi units 
Lowpass Filter Magnitude in dB 
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Figure 7 Using Impulse-invariance 
transformation method and bilinear elliptic 
low-pass filter method to high-pass filter for 
the frequency and amplitude relationship 

From above Example 7 we can find out that 
digital elliptical filter using Impulse- 
invariance transformation method and 
bilinear elliptic low-pass filter method to 
high-pass filter at the confine specifications 



Wp = [0.2 0.55]; Ws = [0.3 0.5] ; Rp = 
ldB; As = 20dB; will produce a small 
ripples beside the band stop regions and 
obtain a clearly slope for the next region, it is 
a good transform results for a band stop 
elliptic filter. 

IV. Results and Discussions 
This article by the above theoretical 
analysis, first by the analog frequency 
conversion from the regularization of the 
classical analog low-pass filter transfer 
function H a (s) start, pre-distorting the cut-off 
frequency, and then apply the bilinear 
transformation equation to produce the 
equivalent digital filter H(z). For the 
frequency conversion of analog design with a 
diameter-degree low-pass cut-off frequency 
in the elliptic filter is normalized elliptic low- 
pass filter Q from the start, and with 
replacement ofs/Q . Using the same method 
can be normalized low-pass filter, converted 
to other types of frequency selective filters, 
such as high-pass, band-pass, and band-stop. 
This series of computer simulation using 
MATLAB software for analysis of related 
functions and a graphical plot of the 
conversion result, as shown in Figure 1-6. 
Finally, there is the discrimination, by (a) the 
same impulse-invariance method (b) Bilinear 
transformation method, Laplace equation can 
be converted to the Z-transform and discrete 
systems using MATLAB software to analyse 
the selected time interval of its calculation 
and the results plotted graphically in Figure 7, 
the resulting two kinds of digital signal 
conversion of a graphic comparison of the 
results are very similar. This calculation by 
the simulation system using MATLAB 
software for analysis of related functions and 
the results plotted graphically, and to get 
good simulation results. Found in Figure 3 of 
the high-pass digital filter of the amplitude of 
the Bilinear transformation by Impulse- 
invariance produce more amplitude than 
higher-up about 25%. 
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V. 



Conclusions 



(1) This analog elliptic filter system, through the establishment of 
model. Analog signals with the proposed digital signal equation, 
obtained using MATLAB analysis and calculation of the order of 
the impulse response of different situations, and once again 
confirmed that the number of higher order than those of their 
turning point clear. 

(2) for discrete processing and digital filter system to get, by (1) 
Pulse constant conversion method (2) Bilinear transformation 
method and other two conversion methods, Laplace equation can be 
converted to use the Z- conversion of discrete digital system, and 
selected using MATLAB software to analyze its time interval is 
calculated and obtained good simulation results, and bilinear 
transformation method to get the same conversion than the pulse 
method has a better conversion results. 
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